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impurity  potential  is  of  the  &-function  type;  i.e.,  it  has  only 
1  matrix  element  between  Wannier  functions,  a  diagonal  one  refer¬ 
ring  to  1  site.  We  show  that  this  result  also  holds  for  a  spheri¬ 
cal  hand  in  3  dimensions.  Our  2nd  result  is  that,  for  the  Coulomb- 
impurity  potential  screened  by  the  static  dielectric  constant,  and 
the  above  spherical  band,  there  are  hydrogen-like  localized  states 
built  out  of  states  near  the  subsidiary  minimum,  even  when  all 
powers  of  k  in  E(k)  are  taken  into  account  in  the  equation  for  the 
envelope  function.  The  deviation  of  the  impurity  potential  from 


slow  variation  causes  a  long  lifetime  for  decay  of  the  localized 


state  into  conduction  states  of  the  same  energy. 

~  -8 

shallow  impurity  state ,  the  lifetime  is  =  10 


For  a  typical 
>9 

10  sec. 
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Recently  the  existence  of  discrete  localized  electronic  states, 
arising  because  of  a  localized  impurity,  and  coincident  in  energy  with 
states  in  an  energy  band  has  been  invoked  in  order  to  explain  a  number 
of  observations  in  semiconductors.^  These  states  are  of  two  slightly 
different  types. 

The  first  type  is  associated  with  the  valence  band  systems  in  silicon 
and  germanium.  Spin-orbit  coupling  causes  a  partial  breakdown  of  what 
would  otherwise  be  a  sixfold  degeneracy  at  jt  =  0  into  two-  and  four-fold 
degeneracies,  the  fourfold  degenerate  band  having  the  higher  energy. 

If  the  perturbation  potential  caused  by  a  localized  impurity  such  as  a 
type  III  or  type  V  substitutional  atom  is  sufficiently  weak  (i.e.  its 
matrix  elements  are  not  comparable  in  size  to  the  spin-orbit  splitting 
at  k  =  0),  and  it  varies  slowly  over  distances  of  unit  cell  dimensions, 
so  that  effective  mass  theory  holds,  one  can  imagine  that  the  two-fold 
and  the  four-fold  degenerate  bands  separately  give  rise  to  hydrogen-like 
localized  states.  In  this  case  a  localized  state  arising  from  the  two¬ 
fold  degenerate  band  will  have  the  same  energy  as  some  continuum  state 
of  the  four-fold  degenerate  band. 

A  second,  similar  situation  arises  when  one  considers  the  conduction 
band  of  germanium,  which  has  four  absolute  minima  at  the  ends  of  the 
[l,l,l,]  axes  in  the  Brillouin  zone,  and  a  subsidiary  minimum  at^ic  =  0. 

As  above,  one  may  suspect  that  localized  states  arise  composed  of  Bloch 
states  with  vectors  near  the  subsidiary  minimum,  in  addition  to  the 
localized  states  that  arise  from  Bloch  states  near  the  absolute  minima. 

The  weak  point  in  such  speculations  is  that  the  usual  effective  mass 
theory,  upon  which  they  are  based  (a  form  of  perturbation  theory)  ignores 
the  existence  of  possible  continuum  states  degenerate  with  the  localized 
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state  that  is  created.  Although  a  slowly  varying  potential  has  very 
small  matrix  elements  between  states  with  widely  different  k  values,  a 
small  or  vanishing  unperturbed  energy  difference  between  these  states 
may  cause  appreciable  mixings  of  these  types  of  unperturbed  stages  in  any 
state  describing  the  perturbed  system.  As  an  example  of  such  mixing  we  may 
mention  that  for  a  strongly  localized  impurity  potential  perturbing  a  mon¬ 
atomic  one -dimensional  solid  having  an  energy  band  with  a  subsidiary  minimum, 
2 

Slater  and  Koster  have  shown  that  no  completely  localized  state  exists 
associated  with  the  subsidiary  minimum  regardless  of  the  strength  of  the 
localized  impurity  potential.  The  purpose  of  this  paper  is  to  show  that  if 
certain  restrictions  of  the  effective  mass  approximation  are  removed,  allow¬ 
ing  one  to  effectively  include  interactions  of  all  Bloch  waves  in  an 
energy  band,  and  if,  in  addition,  the  effective  potential  caused  by  the 
localized  impurity  is  still  slowly  varying,  rigorously  localized  states 
exist  associated  with  a  subsidiary  minimum  of  the  band.  The  remainder 
of  the  perturbation  potential  is  estimated  to  be  responsible  for  a  life¬ 
time  long  compared  to  a  typical  lifetime  caused  by  the  electron  lattice 
interaction. 

Let  us  first  review  some  aspects  of  Slater  and  Koster 's  one -dimensional 
calculation.  They  consider  an  infinite  one -dimensional  monatomic  crystal 
which  contains  a  localized  impurity.  The  wave  function  \[r(x)  for  the 
perturbed  problem  is  expanded  as  a  linear  combination  of  Wennier  func¬ 
tion  0n(x-pa)  for  the  different  bands  defined  by  the  unperturbed,  peri¬ 
odic  Hamiltonian.  Thus,  \jf(x)  =p^n^n(p)0n(x  “  Pc)  where  p  is  an  integer, 
n  is  the  band  index,  and  a  is  the  lattice  spacing.  The  sum  over  p  extends 
from  -oo  to  +oo.  The  simple  case  is  treated  in  which  the  bands  are 
decoupled,  i.e.  the  impurity  potential  energy  has  no  interband  matrix 
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element  between  Wannier  functions ,  and  further,  for  any  one  band  only 

one  matrix  element  of  the  perturbation  is  non-zero:  /0  (x-pa)  V(x)0(x-pa)  = 


^o^qo^po*  cons*^-er  a  simple  energy  band  having  a  subsidiary  minimum 
in  E(.k)  defined  by  the  unperturbed  problem.  The  equations  defining  the 
amplitudes  U(p)  for  the  unperturbed  problem  are  taken  to  be 


o  -(-t  f  t0)  U(f>)  ju(fi  )  i-  U(p  -  o] 

[U(ftt)  +  ao  -  - 


(i) 


-  -i  o  i 
>  '  j  j 
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where  £  =  /0(x) H  (x)0(x-/o)dx.  Only  enough  terms  £  have  been  kept 

u  O  3/ 

so  that  the  subsidiary  minimum  can  occur  in  E(k).  The  solutions  of 
this  problem  are 


and 


~  £ 


+  2? 


r-  A 


+  2  <5,  2 


(2> 


If  f1>0lBd4i2<-  E(!  :)  exhibits  a  subsidiary  minimum  as  shown 

in  Figure  1. 

The  perturbed  problem  has  solutions  that  are  either  symmetric,  U(p)=U(-p) 
or  antisymmetric  U(p)=-U(-p),  in  p  and  because  of  the  simple  nature  of  the 
impurity  energy  V  it  only  affects  the  symmetric  solutions.  The  equations 
for  the  perturbed  problem  are  the  same  as  for  the  unperturbed  problem  except 
for  p  =  0,  +  1.  By  making  use  of  the  remaining  unperturbed  equations 
and  the  boundary  condition  that  the  wave  function  U(p)  must  not  become 
indefinitely  large  as  p  gets  large  we  find  that  we  may  write 
U(p)  =  Ae  ^  +  B  cos(kp-oc)  for  p  >  0,  when  we  are  in  the  energy  range  A 
in  Fig.  1,  in  which  localized  states  associated  with  the  subsidiary  minimum 
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would  be  expected  to  occur.  Here  k  =  and  iy=  k^  are  given^for  E 
fixed,  by  the  two  solution  for  k'  of  the  dispersion  relation  1^..  (2)»  The 
coefficients  A  and  B  are  now  determined  by  the  homogeneous  equations 


UO)  '  u(-t)  ~o  , 

K  O(o)  t  £?.  (  Of?)  -  0(-s)]  =  o 


(3) 


J 


where  in  these  equations  U(-p)reans  U(p)  for  p  >  0  evaluated  formally 
for  negative  values  of  p.  The  existence  of  a  localized  solution  reqiures 
B  =  0,  and  it  is  easy  to  show  that  no  such  solutions  exist  in  the  energy 
range  A. 

Let  us  now  consider  a  three-dimensional  problem  analogous  to 
Slater  and  Koster's  one -dimensional  model.  We  consider  an  infinite 
three-dimensional  crystal  with  one  substitutional  localized  impurity. 

We  assume  that  there  is  an  energy  bend  of  the  pure  infinite  crystal  which 
is  far  from  other  bands  and  has  a  subsidiary  minimum,  and  that  there  are 
no  matrix  elements  of  the  impurity  potential  energy  between  states  from 
different  bands.  The  Schrodinger  equation  in  the  Wannier  representation 
cannot  be  solved  exactly  in  this  case  and  it  is  more  convenient  to  work 
with  the  'bummation  equation"  form  of  the  difference  equations  for  the 
Wannier  function  amplitudes  as  derived  by  Slater  and  Roster, The  summation 
equation  is 


i 


1 

‘A 


G. 


1 


i/  „  64. 


(4) 


Here,  U.  (=  U(H))  is  the  amplitude  of  the  Wannier  function  centered  at 
i 

lattice  point  S.  Vik  =  /0  (r-R^)  V(r)0Cr-Hj.)dJr  is  the  matrix  element 
of  the  impurity  potential  energy  between  Wannier  functions  on  sites  H. 


and 


JR^,  and  Q  is  the  Green’s  function  defined  by  the  unperturbed 
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periodic  potential  problem: 


Q'-Kl 


7 


'  £  ~  E(A)  ~~ 


b.  2. 


RT ryj  '  E-T$J  (5) 


-ik  •  ft 

where  E(k)  =  £  £^e  ^ -nn.  Ws  wall  now  study  a  model  which  approximates 

our  problem  but  for  which  analytic  answers  can  be  obtained.  We  consider 

a  spherical  band  (E(k)  only  where  k  =  |  k  !)  find  also  consider  the  Brillouin 

zone  to  be  spherical.  A  typical  form  for  E  might  be  E(k)  =  £\^cosk  + 

£  cos2k  and  k  =  n  would  be  the  radius  of  the  zone.  We  keep  two 
2  max 

terms  in  E(k)  in  order  that  there  be  a  subsidiary  minimum,  and  will 
evaluate  the  Green's  function  for  an  energy  below  the  subsidiary  mini¬ 
mum,  but  above  the  bottom  of  the  band.  Because  we  are  only  interested 
in  a  qualitative  question,  namely  the  admixture  of  continuum  and  bound 
components  of  states  lying  within  an  unperturbed  energy  band  but  below 

2  4 

a  subsidiary  minimum,  we  may  further  approximate  E(k)  by  E(k)  =  2ctk  -  0k 
where  both  a  and  0  are  greater  than  zero,  thus  ensuring  the  existence 
of  a  subsidiary  minimum.  We  are  particularly  interested  here  in  the 
case  of  the  very  localized  inpurity  potential 

Kj  -  K  &(&,*)  -  <«> 

In  order  to  study  this  case  we  must  have  a  finite  value  for  G  •  In 
*  oo 

the  case  that  no  subsidiary  minimum  exists  and  E(k)  may  be  satisfactorily 

2 

approximated  by  E(k)  =  2ak  it  is  necessary  to  treat  seriously  the 
existence  of  a  high  wave  number  cutoff  (the  Brillouin  zone  boundary) 

in  the  integral  for  G  .  Failure  to  do  so  leads  to  an  infinite  value 

oo 

for  G  .  However,  in  the  case  we  are  treating  such  a  cutoff  is  not 
oo 
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necessary  to  avoid  a  divergence  of  G.  ,  and  for  mathematical  simplicity, 

oo 

ve  extend  the  range  of  the  integration  in  k  space  to  include  all  of 
space.  In  this  vay  we  find 


2  2  2  4 

where  b  and  -c  are  solutions  of  E  -  23k  +  |3k  =0.  We  nay  of  course 

add  to  this  result  any  solution  of  the  homogeneous  equation,  i.e.  the 
Schrodinger  equation  in  the  Wannier  function  representation  for  the  peri¬ 
odic  potential  problem  with  energ>  E.  We  may  usefully  consider  such  solu- 
tions  to  be  e  J  where  E(k)  =  E.  We  thus  consider  Green's  func¬ 

tions  of  the  form 


r  ° 

<cr. 


'i 


I 


which  give  rise  to  solutions  of  a  perturbed  problem  resembling  as  closely 
as  possible  unperturbed  Bloch  states.  If  the  impurity  potential  is  of 
the  fonn  (6)  then  the  equation  (5)  gives  rise  to  the  eigenvalue  condition 


or 


c 

‘f-rp(^c^) 


+  K(k)- 


(9) 


We  see  that  for  any  energy  that  would  give  rise  to  a  discrete  level 
associated  with  the  subsidiary  minimum  we  may  choose  K(lt)  to  satisfy 
equation  (9).  Thus  we  see  that  no  completely  bound  state  can  exist  for 
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a  very  localized  potential.  Even  if  equation  (9)  could  be  solved  with 
K(k)  =  0,  the  resultant  solution  (for  a  single  energy)  would  contain 
almost  equal  portions  of  localized  and  continuum  parts.  We  have  thus 
extended  Slater  and  Koster's  result  for  the  absence  of  localized  levels 
associated  with  a  subsidiary  minimum  of  an  energy  bend  and  a  strongly 
localized  impurity  potential  to  a  three  dimensional  case. 

Let  us  now  for  comparison  examine  the  case  of  a  Coulomb  impurity 
potential  screened  by  the  static  dielectric  constant.  The  main  point 
of  our  work  is  to  show  that  this  potential  gives  rise^for  a  spherical 
band  rigorously.  in  the  limit  that  the  potential  is  slowly  varying,  and 
approximately  in  other  cases,  to  bound  states  which  are  coincident  in 
energy  with  continuum  states  of  the  perturbed  problem.  For  this  purpose 
it  is  convenient  to  write  the  Schrodinger  equation  in  the  Bloch 
representation: 

Z  ('»&/  H0+Ul'»'£)  ^  -  z  Afe) 


9l 


(10) 


2 

Here  U(r)  =  e/£r(6is  the  static  dielectric  constant)  and 

Hq  =  V(r)  where  V(r)  is  the  periodic  potential.  U(r)  is, 

of  course,  an  approximate  impurity  potential  energy  which  has  been  justi- 

3 

fied  for  semiconductors  by  many  workers  using  recent  many-body  theory. 

4 

Following  Kohn  we  now  write  these  equations  in  a  form  which  separates 
the  isolated  band  (subscript  0)  under  discussion  from  all  other  bands  ' 


O 
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+  2 

■-n'i-O 

fj 


(°&iuh'£)  -° 


V  ) 
/ 


(11) 


2.  h&luloA‘)AM 

+  &)*«&)•  °  *°> 

fl 

where  A (k )  is  the  amplitude  of  the  Bloch  function  with  wave  vector  k^ 
in  hand  n.  Here  we  have  used  the  fact  that 


M  /  K>  h'*' )  =  ,>  £ 


'ji^1  ■ 


The  matrix  elements  of  the  perturbation  arc  given  by 


\ 

ia  1) 


where  k  is  a  reciprocal  lattice  vector  and  Cv  ,  .  is  the  vth  coef- 

v  nkjn'&' 

ficient  in  the  fourier  expansion  of  the  product  of  the  periodic  parts 

of  the  Bloch  functions  \lr  (r)  and  \k  ,  (r).  In  particular  C°  ,  =5 

nk  Tn'kf  *  nk:n’k  nn’ 

V*  vt*.  •*», 

The  Coulomb  potential  is  slowly  varying  except  in  the  cell  containing 
the  impurity  atom.  In  the  familiar  situation  of  a  band  with  no  subsidi- 
ary  minimum  it  is  assumed  that  only  Bloch  functions  with  wave  vectors 
near  the  minimum  will  contribute  appreciably  to  any  bound  impurity  state 
wave  functions.  Thus  all  terms  in  Eq.  (12)  withj£v  =  0  are  neglected, 
and  all  interband  terns  are  at  first  neglected  because  the  impurity  poten¬ 
tial  matrix  elements  are  much  smaller  than  typical  interband  energies.  Further- 


*  CL  ,  2 

more  iQ£)  is  approximated  by  E^  +  ^  (k  ~  k^)  , 


an  expression  valid  near 


G 


the  band,  minimum  A,  The  restriction  that  k  lie  within  the  Brillouin  zone 
is  then  lifted  and  the  usual  slowly  varying  hydrogen  atom  type  envelope 
functions  result^  together  with  a  hydrogen  at. ora  energy  spectrum,  for  the 
bound  levels.  Finally  one  can  check  that  the  envelope  functions  only 
involve  small  k.  values,  thus  establishing  the  consistency  of  the  solution. 

When  a  subsidiary  minimum  exists  it  is  not  clear  that  all  K  ^  0 
can  be  neglected.  Let  us  nevertheless  make  the  same  assumptions  for  the 
subsidiary  minimum  problem  as  were  made  in  the  case  in  which  there  is 
only  an  absolute  minimum  and  examine  the  justification  later.  We  neglect 
all  and  interband  effects,  and  we  extend  the  range  of  k^  out  of 

the  Brillouin  zone  into  all  reciprocal  space.  However  we  now  use  an 
E(fe.)  that  allows  a  subsidiary  minimum.  In  addition  we  assune  a  spheri¬ 
cal  band  in  order  that  we  can  solve  the  equations  for  the  envelope  func¬ 
tion.  The  impurity  equation  in  space  then  becomes 


(Ra-eKWdi 


c 


'i  2 


Aik  )d£  ~0)  (13) 


where  the  sum  over  kr  in  (ll)  has  been  changed  to  an  integral. 

We  take  E(k)  =  +  £”.  cosk  +  £  cos 2k  where  the  conditions 

'  o  1  2 

<T,  +  4£  <  0  and  £.  >  0  establish  a  subsidiary  minimum  at  k  =  0  and 
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an  absolute  minimum  at  k.  =  n.  Here  k  is  the  magnitude  of  k. 

If  we  now  return  the  problem  to  "real"  space  by  multiplying  by 
and  integrating  over  reciprocal  space,  we  find,  the  following 
equation  for  the  envelope  function  x(r^) 

+ s  vrf-  n/rjxfch  -a$h) 


-li- 


where  §  =  f)  =  e2/^,  and  £  =  (E-  E^IS  x*  Here  we  understand 

the  cos  terms  to  mean  their  power  series*  Because  the  Hamiltonian  is 
spherically  synmetric  we  can  write 


We  then  find  that 


Let  us  first  treat  the  case  l  •  0.  We  let  R(r)  =  F(r)/r  and  find  that 
the  equation  tor  F  becomes 


<50 


z 


F 


•2-h 


*0  M! 


_  -  r 

(?■*)!  ri  > 


+  $  Z  S¥ 


(15) 


e 


where  F 


2n 


d^V 


dr 


2n 


This  equation  is  a  generalization  of  the  confluent  hypergeometric 
equation5,  and  it  is  thus  not  surprising  that  when  we  try  a  solution  of 
the  form  R(r)  =  e  the  resulting  equations  for  a. 


OsCtgJL  -f  £  &  j  (a) 

—  JL*  ~  "7  (y>\ 
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have  a  solution  with  a  real.  These  equations  arise  when  the  coefficients 
of  the  different  powers  of  r  are  set  equal  to  zero.  Eq.  l6  (a)  gives 
the  energy  when  a  is  known  and  Eq.  1 6  (b)  gives  a  in  terms  of  the  para¬ 
meters  of  the  problem,  including  the  strength  of  the  impurity  potential. 
Equation  l6  (b)  has  two  solutions  for  a,  one  real  and  positive,  and  the 
other  of  the  form  a  -  in  +  X,  where  y  is  positive.  These  solutions  cor¬ 
respond  to  bound  states  splitting  off  from  the  band  at  k  =  0  and  k  =  n 

AVt 

respectively. 

For  a  or  x  small  the  equations  may  be  solved  approximately  yielding 
the  usual  effective  mass  results :  _ 


o<  - 


=  -2L~.  e.  =  -r-/  +  ZL 


i  +*  £  J 


V  =  -_JL 


-I  ++S 


)  £y  -  -S+l 


*|l  -LV§I  J 

2  (17) 

-+  2L 


It  can  be  seen  that  for  the  form  of  E(k)  that  we  used  a  >  x  no  matter 

what  the  values  of  q  and  5  but  that  equality  of  these  quantities  is 

approached  when  5  becomes  large.  Bound  states  of  higher  energy  for 

2  =  0  can  undoubtedly  be  constructed.  For  2  /flit  is  not  possible  to 

solve  the  differential  equation  exactly.  However,  for  large  r  all  con- 

tributions  from  2{2+±)/r  in  Equation  l4  are  negligible  so  that  the 

■ocr 

asymptotic  solution  e  ,  where  a  has  a  positive  real  part  is  again 
possible. 

Our  result  shows  that,  to  the  extent  that  the  potential  energy  in 
the  equation  for  the  envelope  function  can  be  assumed  to  be  Coulomb-like j 
bound  impurity  states  of  the  effective  mass  type  may  be  associated  with  a 
subsidiary  minimum  in  a  spherical  band.  In  contrast  with  the  case  of 
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the  very  localized  impurity,  it  is  not  true  that  all  states  with  energies 
in  the  unperturbed  band  region  contain  admixtures  of  continuum  states. 

It  is  reasonable  to  expect  that  this  result  does  not  depend  critically 
on  the  use  of  a  spherical  band  but  we  have  not  studied  this  point  in 
detail. 


Of  great  interest  is  the  possible  effect  of  a  deviation  of  the  poten¬ 
tial  that  should  be  inserted  in  the  envelope  equation  from  the  Coulomb 
form.  In  this  connection  we  remember  that  there  are  two  effects  to  con¬ 
sider:  the  neglect  of  jj  f  C  and  the  extension  of  the  Brillouin  zone 
to  all  of  space.  These  effects  are  indeed  connected.  Consider  first 

the  effect  of  the  extension  of  the  band  to  all  reciprocal  space  on  the 
2 

term  1/|  £  -  1 1  ■  If  the  band  were  not  extended  we  would  be  omitting 
high  k  -  k',  thus  putting  a  bottom  on  the  potential  well.  This  should 


not  matter  as  long  as  no  bound  states 

Now  consider  the  term  proportional  to 

of  this  term,  CV  •  „  is  not  likely 
OK*  OK 


lie  near  the  bottom  of  the  well. 

p 

l/J'k  -  k*  -  K^l  .  The  coefficient 
to  be  greater  than  unity  but  can 


be  counted  on  to  be  of  order  unity  until  v  gets  large,  in  which  case  no 


large  contributions  are  expected  anyway.  We  will  neglect  the  k  and  k* 


dependences  of  Cv  and  set  this  coefficient  equal  to  unity.  We  are 

OK  •  OK 

thus  led  to  terms  in  the  potential  energy  of  the  form 


-14- 


If  we  extend  the  boundary  of  the  Brillouin  zone  to  infinity  this  terra 


gives 


e 


r 


d 3 (a  -a  -  ^ ; 

J  ( - 

L  K  -n. 


eT  &  ~  ^  * 

^  v_. 


/Xs 


Including  both  ft  and  -ft  ,  ve  get  a  contribution  to  the  potential  of  the 
v  -**v 


form 


V  <=<  c^i  iO  ’& 

n 


(19) 


Terras  such  as  Eq.  (.1$)  are  smaller  then  the  l/r  we  have  been  using. 

If  the  zone  is  not  extended  to  infinity,  the  integrals  become  very  hard 
to  do,  and  it  is  questionable  whether  the  approximations  of  the  problem 
warrant  such  a  calculation.  This  discussion  shows  that  it  is  not  easy 
to  find  the  proper  potential  to  use  in  the  envelope  function  equation 
when  the  real  potential  (for  low  lying  states)  in  the  Scbrodinger  equa¬ 
tion  is  the  Coulomb  potential.  However,  it  is  also  indicated  that  the 
Coulomb  potential  is  surely  of  the  right  order  of  magnitude  and  it  is 
a  bit  surprising  that  bound  states  coincident  in  energy  with  continuum 
states  exist  for  this  potential.  Perhaps  there  is  a  mathematical  con¬ 
nection  with  the  fact  that  the  Klein-Gordon  equation,  in  which  the 
kinetic  energy  can  also  be  thought  of  as  a  series  in  powers  of  the 
Laplacian  operator,  has  bound  solutions  for  the  hydrogen  atom. 

A  change  in  the  potential  from  the  Coulomb  form  will  probably  change 
the  wave  functions  enough  to  cause  bound  states  coincident  in  energy  with 
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continuum  states  to  no  longer  be  eigenstates  of  the  impurity  problem, 

but  .if  so  the  lifetimes  of  these  states  will  -probably  be  long.  An 
estimate  was  made  of  the  lifetime  of  a  bound  state  of  the  form  u  (r)e 
with  respect  to  decay  into  continuum  states  of  tne  same  energy,  using 
a  potential  of  the  form  (19)  summed  over  the  smallest  reciprocal  lattice 
vectors  and  a  spherical  band  of  the  form  used  above.  The  result  for 
the  width  £E  of  the  bound  state  is 


AE  Ce  3 m^sV%e)(^)(=<a 

)f ^ 

/ 

1  \zer  - 

he 

'>(20) 


where  B  and  (8  arise  from  expansions  of  products  of  the  periodic  parts 

of  two  Bloch  functions  from  the  some  band  in  terms  of  a  Fourier  series 

6 

in  reciprocal  lattice  vectors,  and  are  of  order  of  unity;  a  is  the 
Bohr  radius,  and  all  other  symbols  have  been  defined  earlier.  Here 
e  a/2  6  is  the  energy  of  a  shallow  impurity  and  we  see  that  £E  is 
smaller  than  this  by  two  factors,  the  first  being  the  ratio  of  the 
impurity  energy  to  an  energy  of  the  order  of  a  band  width,  and  the 

second  the  ratio  of  the  Bohr  radius  to  the  impurity  orbit  radius. 

-8  -9 

Eq.  (20)  gives  a  lifetime  for  a  typical  shallow  impurity  of  10  — 10 

sec. 

In  summary,  we  may  say  that  impurity  levels  with  very  long 
lifetimes  are  associated  with  subsidiary  minima  for  a  Coulomb  impurity 
potential,  although  no  such  levels  exist  for  a  very  localized  impurity 
potential.  It  would  be  interesting  to  know  how  slowly  the  impurity 
potential  must  vary  in  order  that  localized  impurity  levels  be  associated 
with  subsidiary  minima. 
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E(k)  for  an  energy  band  having  a  subsidiary  minimum 


